TWO-DIMENSIONAL INCOMPRESSIBLE VISCOUS FLOW 
AROUND A SMALL OBSTACLE 



D. IFTIMIE, M. C. LOPES FILHO AND H. J. NUSSENZVEIG LOPES 

Abstract. In this work we study the asymptotic behavior of viscous incom- 
pressible 2D flow in the exterior of a small material obstacle. We fix the 
initial vorticity loo and the circulation 7 of the initial flow around the obstacle. 
We prove that, if 7 is sufficiently small, the limit flow satisfies the full-plane 
Navier-Stokes system, with initial vorticity loo + 7<5> where <5 is the standard 
Dirac measure. The result should be contrasted with the corresponding invis- 
cid result obtained by the authors in 1151 . where the effect of the small obstacle 
appears in the coefficients of the PDE and not only on the initial data. The 
main ingredients of the proof are L p — L q estimates for the Stokes operator in 
an exterior domain, a priori estimates inspired on Kato's fixed point method, 
energy estimates, renormalization and interpolation. 



Contents 



Introduction 



^gtirn^tcs_jor_th£_ ^tokc^ semigroup 



^h^^vam^jccrr^ojostack 



^M^^^f^^^^nrjtotics 



5. The impulsively stopped rotating cylinder 



6. Initial-laver and the nonlinear evolution! 



J_^_ G_j£bal^ri-yrncjKmlincar cvolutio 



„ 



Velocity estimates' 



9. Compactness in s pace-time 



10. Passin g to the limit 

11. Conclusion J ~ 
References! 



2 
3 
7 
11 
13 
17 
20 
22 
24 
29 
3G 
38 



2 



IFTIMIE, LOPES FILHO AND NUSSENZVEIG LOPES 



1. Introduction 

The purpose of this work is to study the influence of a material obstacle on the 
behavior of two-dimensional incompressible viscous flows when the size of the ob- 
stacle is small compared to that of a reference spatial scale. More precisely, we fix 
both an initial vorticity uiq, smooth and compactly supported, and the circulation 
7 of the initial velocity around the boundary of the obstacle, while homothetically 
contracting the obstacle to a point P outside the support of wo- The initial vortic- 
ity luq and the circulation 7 uniquely determine a family of divergence-free initial 
velocities Uq with curl Uq = luq and Uq(x) — > at infinity. The size of the support 
of the initial vorticity ojq can be used as reference spatial scale. Let u £ — u £ (x, t) 
be a solution of the Navier-Stokes equations with initial data Uq and no-slip data 
at the boundary of the small obstacle. Our problem is to determine the asymptotic 
behavior of u £ as e — ► 0. We will show that u £ converges to a solution of the 
Navier-Stokes equations in the full plane with initial vorticity luq + "fS(x — P), as 
long as 7 is sufficiently small. 

There is a sharp contrast between the behavior of ideal and viscous flows around 
a small obstacle. In j!5j . the authors studied the vanishing obstacle problem for 
incompressible, ideal, two-dimensional flow. The ideal flow assumption is physically 
incorrect in the presence of material boundaries, and part of the motivation of the 
present work, together with |15j . is to explore more precisely the extent to which the 
ideal flow assumption misrepresents the physical flow in the presence of material 
boundaries. The main result in ^H] is that the limit vorticity in the ideal case 
satisfies a modified vorticity equation of the form ui t + u ■ Vcj = 0, with div u = 
and curl u = uj + jS(x — P). In the ideal case, the correction due to the obstacle 
appears as time-independent additional convection centered at P, whereas in the 
viscous case, the correction appears only on the initial data and gets convected and 
diffused. 

The small obstacle limit is an instance of the general problem of PDE on singu- 
larly perturbed domains. There is a large literature on such problems, specially in 
the elliptic case, see for a broad overview. Asymptotic behavior of fluid flow on 
singularly perturbed domains is a natural subject for analytical investigation which 
is virtually unexplored. The present work, together with ^S], may be regarded as 
a first attempt to address this class of problems. 

There is a natural connection between the approximation problem as we have 
formulated it and the issue of uniqueness for the limit problem. In fact, from a 
technical point of view, our work is closely related to the classical uniqueness result 
due to Y. Giga, T. Miyakawa and H. Osada, on solutions of the incompressible 
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2D Navier-Stokes equations with measures as initial data, see ^3). Some of the 
more striking similarities are: the difficulties with locally infinite kinetic energy, 
the use of L p estimates for the linearized problem and the use of Kato-type norms 
to estimate the nonlinearity. The smallness condition on the mass of the point 
vortices in the initial data, required in the uniqueness result, is closely related to 
our smallness condition on the circulation. 

The remainder of this work is organized in ten sections. In Section 2 we summa- 
rize LP estimates for the time-dependent Stokes problem on exterior domains. In 
Section 3 we formulate precisely the problem we wish to discuss and write uniform 
estimates for the initial data. In Section 4 we study the asymptotic behavior of 
the initial data. In Section 5 we discuss physical motivation for our problem and 
we establish the small obstacle asymtotics for circularly symmetric flows, a linear 
version of our problem. In Section 6 we derive a priori estimates in the initial layer 
for the nonlinear correction term. In Section 7 we deduce global-in-time energy 
estimates for the nonlinear correction term. In Section 8 we put together the esti- 
mates for the linear part with the estimates for the nonlinear correction, obtaining 
a complete set of a priori estimates for velocity. In Section 9 we prove compactness 
space-time, in Section 10 we perform the passage to the limit and in Section 11 we 
add comments and concluding remarks. 

We conclude this introduction with a few remarks regarding notation. Given a 
vector z — (zi,z 2 ) € P 2 we denote its orthogonal vector by z 1 = (—22, Z\). We use 
the subscript c in function spaces to denote compact support, as in , and we 
use standard notation for Sobolev spaces, W k,p , where 1 < p < 00 and k E Z, with 
H k standing for the case p = 2. We use the subscript loc in function spaces X to 
denote functions which are locally in X. In particular, £f oc ([0, 00); W k,q ) denotes 
functions / = f(t, x) S V([0, M]; W k ^ q ) for any M > 0, whereas Lf oc ((0, 00); W k ' q ) 
denotes functions / = f(t, x) e L p ([5, M];W k ' q ) for any 5 > and any M > 0, but 
not necessarily for 5 = 0. Finally, L 2 '°° denotes the Lorentz space of functions / 
whose distribution function satisfies A/ = A/(s) = |{|/| > s}\ = (D(s~ 2 ). 

2. Estimates for the Stokes semigroup 

In this section we will put together several results on estimates for the Stokes 
semigroup on exterior domains. Let us begin by introducing some basic notation. 

Let be a bounded, open, simply connected subset of K 2 with boundary T, a 
smooth Jordan curve. We denote by II the unbounded connected component of 
M 2 \ r. Fix v > and let P denote the Leray projector onto divergence-free vector 
fields on II. Let A = —PA be the Stokes operator on II and denote the Stokes 
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semigroup by S„(t) = e^ tA . Given v € C c °°(n), let v(t,x) = S v (t)v be the 
unique solution of the system 



(2.1) 



dtV — i/Aii = — Vp, in (0, oo) x II 

div v = 0, in [0, oo) x II 

v = 0, on (0, oo) x T 

limi^i^oo v(t, x) = 0, for all t > 

v(0, x) = Vq(x), on {t = 0} x II. 



We denote by X P (H) the closure of the space of divergence- free, C£°(II) vector 
fields with respect to the L p -norm. The Stokes operator in X p generates an analytic 
semigroup of class C° on X P (U), for any 1 < p < oo, see [E], so that, in particular, 
problem 1)2. 1 1> is well-posed in X P (H). 

We will require two kinds of estimates on the Stokes semigroup, L p estimates 
and renormalized energy estimates. We first state the L p estimates. 

Theorem 1. Let 1 < q < oo. Consider v G X q (Tl) and F G L q (Tl; M 2x2 (M)). 
T/ien we /iaue f/ie following estimates. 

(51) Let q < p < oo. There exists K\ = Ki(JI,p, q) > suc/i t/iat 

||5 v (*)«o||t» < Ki(ut)p-*\\vo\\L», 

for all t > 0. 

(52) Let g < p < 2. There exists K% = K2(H,p, q) > such that 

||V^(t)«o||£p < K 2 (ut)-*H-l\\v \\ L „ 

for all t > 0. 

(53) Assume q > 2 and let q < p < oo. Then there exists K% — K%(Jl,p, q) > 



||S„(t)P diw F|| L , < K 3 (vt)-^ + p-i\\F\ 



Li, 



for all t > 0, wi</i i/ie divergence taken along rows of the matrix F . 

This theorem summarizes several results already contained in the literature, 
which we have collated above for convenience. 

Proof. Estimates (SI) and (S2) were proved in [S] [20! (see also |S] for the case 
p = oo). Estimate (S3) follows from (S2) by duality. Indeed, the adjoint of S v (t) 
on X p is again S v (t), defined on X p , with l/p+l/p' = 1 and therefore the adjoint of 
VS v (t) is S v (t) P div . The dependence on the viscosity follows directly by rescaling 
time, since S v (t) = Si(vt) . □ 
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Next we address a renormalized energy estimate for the Stokes semigroup. Our 
concerns include infinite energy solutions to the Navier-Stokes equations whose 
behavior at infinity is 0(1/ |x|). In the following result we will prove that solutions 
to the Stokes system retain the behavior at infinity of their initial data. 

Proposition 2. Let vq be a smooth divergence-free vector field on H vanishing 
at the boundary T. We assume also that vq £ X P (TV) for some p > 2 and that 
Wv Q E L 2 (U). Then 

S v (t)vo -v e C°([0, oo); i 2 (n)) n LL([0, oo); if 1 (II)). 

Moreover the following inequality holds 

(2.2) \\S„(t)v Q - v \\ 2 L2 + v f \\V[S v (t)v - v ]\\l 2 dT < vt\\Vv \\l2. 

Jo 

Proof. Let W = S„(t)v — v . Then W satisfies: 

d t W - vAW = ~Vp + vAv , in (0, oo) x II 

div W = 0, in [0, oo) x II 

(2.3) { W = 0, on(0,oo)xr 
lim^i^oo W(t, x) = 0, for all t > 
W{0,x)=0, on{<^0}xn. 

It is well-known that 1)2.. admits a unique solution W in C°([0, oo); L 2 (II)) (~l 
^ oc ([0,oo);i/ 1 (n)), see, for instance, Theorem III.l.l in The fact that W - 
W = follows from the well-posedness of l|2.1|l in X p . The standard energy estimate 
gives (ESI. □ 

One consequence of the nontrivial topology of II is the existence of harmonic 
vector fields, i.e. divergence-free and curl-free vector fields which are tangent to T 
and vanish at infinity. We denote by i?n the unique harmonic vector field on the 
exterior domain II which satisfies the condition 

Hu ■ ds = 1, 

r 

where the contour integral is taken in the counterclockwise sense. It is an elementary 
application of Hodge theory that the vector space of these harmonic vector fields 
on II is one dimensional, and we can use Hu as a basis. In the case where II is the 
exterior of the unit disk centered at the origin, we will denote Hu simply by H , 
and we have: 
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We will require detailed information on the behavior of Hu both at infinity and 
near the boundary T, which we obtain by means of a conformal mapping. We 
denote 

U = {\x\>l} 

and switch to complex variables notation in the result below. 

Lemma 3. There exists a smooth biholomorphism T : II — * hi, extending smoothly 
up to the boundary, mapping T to {\z\ — 1}. Furthermore, there exists a nonzero 
real number (3 and a bounded holomorphic function h : IT — > C such that: 

(2.5) T(z) = /3z + h(z). 
Additionally, 

(2.6) h'(z)=o(jLY m |*| -oo. 

This Lemma is an excerpt from |15|. Its proof is an exercise in complex analysis. 
It was observed in JHj (see identity (2.10) in JH]) that 

1 DT t (x)(T(x)) x 

From Lemma |3 we see that \Hn\ is C(l/|x|) for large \x\. This implies that Hu 
belongs to the Lorentz space L 2 '°°(n). 

We close this section with an estimate for the Stokes semigroup acting on infinite 
energy initial data. 

Proposition 4. Let 2 < p < oo and let v £ L 2 '°°(n) n X P (U). There exists a 
constant K$ > such that 

\\S u (t)v \\ L * < Ksivt^-iWvoWv^. 

In particular, this estimate holds true for vq = Hu (x) . 

Proof. This estimate is contained in Proposition 2.2, item (4), of ^7]. To see that 
it holds for Hu, we first show that Hu £ X P (H) for any p > 2. This is easy to 
prove in the case IT = U because, for any function ip £ C£°((0, oo)), ip(\x\)H(x) is 
smooth, compactly supported and divergence-free, and, by taking ip £ a sequence 
of cutoffs for the interval (1 + e,l/e), it is easy to see that (p £ (\x\)H — > H in 
L p , for p > 2. For general II, we use the conformal mapping T, approximating 
Hu by ip £ (\T(x)\)Hu(x), where <p e is the same family of cutoffs used in the case 
of the exterior of the disk. This strategy works because (p e (\T(x)\)Hu(x) is also 
divergence-free. □ 
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3. The evanescent obstacle 

The purpose of this section is to set down a precise statement of the small 
obstacle problem. Many of the key issues regarding the small obstacle limit and 
incompressible flow have been discussed in detail in ^H], so that we will focus on 
issues specifically related with viscous flow and briefly outline the rest. 

As in ^I)], nx w o G C^°(1R 2 ) and assume that the origin does not belong to the 
support of luq. Let il be a bounded, open, connected and simply-connected subset 
of the plane whose boundary L is a C°° Jordan curve. The evanescent obstacle 
is the family of domains sil, with < e < eq. The parameter Sq is chosen small 
enough so that the support of u>o does not intercept efl for any < e < £q. 

Fix < e < Eq. Let LT £ = M 2 \ e£l and T e = dH e . We use the conformal mapping 
T : IT U, given in Lemma to define a family of smooth biholomorphisms 

(3.1) T e = T e (x) = T(~y 

Throughout we write H £ for Hjj e and G £ = G £ (x,y) will be the Green's function 
of the Laplacian in II e . Let K £ (x,y) — V^G £ (x,y) be the kernel of the Biot- 
Savart law on II e and denote the associated integral operator by / i— > = 
J n K £ (x,y)f(y)dy. Both K £ and H £ are related to Ku and Hu respectively, 
through the conformal mapping T £ , in a way which was made explicit in [15] . The 
relevant fact is the way that both the Biot-Savart kernel and the basic harmonic 
vector field scale with e, see identities (3.5) and (3.6) in |15| . 
Fix a G M and let 

(3.2) u £ Q = K £ [u ] + aH £ . 
We consider the problem 

d t u £ + u £ ■ Vu £ - vAu £ = -Vp £ , in (0, oo) x U E 

div u £ = 0, in [0, oo) x II e 

(3.3) \ u £ = 0, on(0,oo)xL e 
limui—joo u £ (t, x) = 0, for all i > 
w e (0,x) = u%{x), on {t = 0}xn £ . 

We begin by observing that G L 2 '°°(II e ) nL p (II e ) for any 2 < p < oo. Indeed, 
Uq is smooth, and therefore locally bounded, so that we only require knowledge on 
the behavior of Uq at infinity. By Lemma and identity (|2.7|l \H £ \ has 0(1/ |x|) 
behavior as |a;| — > oo, and therefore it belongs to L 2 '°°(n £ ) n L p (Tl e ) for any 2 < 
p < oo. In fact, the L 2 '°° bound on \H £ \ is independent of e as can be readily seen 
by rescaling to a fixed domain and using that Hu belongs to L 2 co . In ^Hj it was 
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shown that |fT £ [o;o]| has behavior 0(l/|a;| 2 ) at infinity (see estimate (2.8) in [15j ) 
and therefore it belongs to L P (H £ ), for any p > 2, and, in particular, to L '°°(II e ). 

Global-in-time well-posedness for problem l|3.3|l was established by Kozono and 
Yamazaki in |17|. The existence part of Kozono and Yamazaki's result requires 
that the initial velocity satisfy a smallness condition of the form 

limsupi? | {a; G II £ | |ttg(»)| > R}\ 1/2 < 1. 

Since Uq is bounded, the limsup above is always zero, for any e > 0. Uniqueness 
holds for divergence- free initial data in L 2, °°+X p without any additional conditions. 

The evanescent obstacle problem consists of understanding the asymptotic be- 
havior of Kozono and Yamazaki's solution u £ (x,t) for small e. More precisely we 
will show that, under appropriate assumptions, u £ has a limit, and we will identify 
an equation satisfied by this limit. 

Fix <p : K — ► [0, 1] a smooth, monotone function such that ip(s) = if s < 2 and 
(p(s) = 1 if s > 3. For each e > and A > we introduce the adapted cut-off 
functions: 

(3.4) ^(ijEpgrei), 

Note that the cutoff function ip £,x vanishes in a ball of radius 0(A) and it is identi- 
cally equal to 1 outside a larger ball of radius 0(A), for large A. Furthermore, the 
radii of the annulus where tp £ ' X is not constant can be made independent of e. This 
follows easily from the fact that T is asymptotically affine at infinity, see 12. 5[l . 

We will now introduce a pair of parameters that are useful to describe the as- 
ymptotic behavior of Uq when e — > 0. Consider 

(3.5) m = J uJodx and 7 = j) u s Q ■ ds. 

By Stokes' Theorem we have that 7 = a — m, and therefore, the circulation 7 does 
not the depend on e, see the proof of Lemma 3.1 in |15| . 

For each A > 0, we introduce a convenient decomposition of the initial velocity 

as 

= b o + io + °o. 

with 

=K £ [u ] +m(l-<p s > x )H s , 

it EE 7 (1 - ^)H £ , 

and 

= aif £ ' X H £ . 
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We need to understand the behavior of each of the components of this decom- 
position, in the limit e — > 0. This is the content of our next result. The proof uses 
a large part of the work done in JS] . 

Lemma 5. There exists \q > 0, independent of e, for which 1 1 1>5 1 1 ? l|iolU P2 
and ||og||iP3 are uniformly bounded in e, for any 1 < p± < oo, 1 < P2 < 2 and 
2 < P3 < oo. The vector fields bg, ig and Og are divergence-free, the first two are 
tangent to T e and the last one vanishes on T e . Moreover, || VogH^m^ is bounded 
independent of e and 



(3.6) 



— > as e — > 0, 

i 2 (n e ) 



An 

where H = x^ / (2ir\x\ 2 ) and (3 is as in Lemma\Q We also have that 

WUl^) <C|log £ |i 

Proof. Choose Ao such that the radii of the annulus where (p £ - x ° is not constant are 
uniform in e < e . 

The L°° bound on bg comes from Theorem 4.1 in U^J- The L Pl bound on bg, 
1 < pi < oo follows from the local L°° bound above, together with two facts: (1) 
m(l — ip £ ' X °)H e has support in a compact set independent of e and (2) |A" e [o;o]| = 
0(l/|:r| 2 ) at infinity, uniformly in e. As mentioned previously, fact (2) is estimate 
(2.8) in P2]. 

The L P3 bound on Oq follows from the fact that ip £,x ° is constant outside an 
annulus independent of e, from formula i|2. 7fl . from the scaling H e (x) — ^Hn{x/e) 
and from the behavior of T far from the obstacle given by Lemma |3 

For ig, both the logarithmic estimate and the L P2 estimate follow from adapting 
the argument used for estimate (3.7) of JHj m a straightforward manner. 

To estimate on Vog we observe that |Vog| = C(l/|x| 2 ) near infinity, uniformly 
in e. This estimate easily reduces to an estimate on DH e , which in turn reduces 
to calculating derivatives of the conformal mapping T using H2.5|l . 

Finally, (|3.6|l reduces to showing that H — H £ goes to zero in L 2 near infinity, 
which can be shown by a computation similar to the one carried out in the proof 
of Lemma 4.2 in |15j . 

□ 

In the remainder of this article, we will fix Ag, independent of e, as in Lemma 
thereby fixing the bounded, inner and outer parts of the initial velocity, bg, ig and 

Let us denote the Stokes semigroup on n e by S £ (t), so that <S^(i)[vp] is the 
solution to the Stokes system H2.lt on II e with initial data Vq. We introduce the 
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notation t £ = t £ (x) = ex, the contraction by e. We observe the following funda- 
mental relation between the Stokes system on LT e and on II: 

(3.7) (Sl(tM o t £ ])(x) = (St(sHM])(sx), if x G II. 

Our strategy to study the small obstacle limit begins by considering the solu- 
tion u £ of (|3.3[1 as a perturbation of V s = S £ (t)uf). The first thing we require is 
information on v £ which we deduce in the result below. 

Lemma 6. Let h £ = Sl{t)h%, i £ = S^(t)% o £ = S £ (t)o £ a and let 2 < p < oo. Then 
there exists a constant K — K(p,luq) > such that for any e > we have: 

(i) l|b £ || LP( n E ) <K{yt)*-$, 

(ii) l|i £ ||L P( n e ) < K\»i\iyt)>-* , 

(iii) ||o e || LP( n e) < K\a\{vt)»-i. 

Proof. By (JSTJ) we have that h £ {e 2 t,ex) = (Sl(t)[b e o t £ ]){x), for x € n. Now, by 
Theorem^ item (SI), it follows that there exists K\ > such that 

ll#(t)[b§ ° r e ]|U, ( n) < K^-^m o r £ || i2( n). 
Item (i) above follows from this estimate, together with i|3.7|) and the fact that 

im°^ £ ||L2(n) = J||bS||i,2 ( n c ) < C^, 

where we have used Lemma [5] in the last inequality. Items (ii) and (iii) follow in an 

analogous manner using Proposition ^together with the fact that 

1 £ | 7 | 
II "-o T E |U 2 -~(n) = -||iollt 2 '~(n e ) < C— 

and 

||ogor £ || L 2,oo ( n) = i||o^|| L 2. TO( n e) < 
We have used the scaling H e (x) = (l/e)H 1 (x/s) above, see identity (3.6) in |J2]. □ 

Remark 7. Using the rescaling (|3.7|) we may deduce that the estimates (SI), (S2) 
and (S3) in Theorem^are valid in II e with constants K\, K% and K3 independent 
of e. 

We will conclude this section with an observation on the amount of vorticity 
generated at the boundary in the initial layer. This is a "fixed e" calculation, 
before we take the vanishing obstacle limit. We denote vorticity associated to the 
velocity u £ , at time t, by oj £ — uj £ (t, ■) = curl u £ {t, •). Let us recall the discussion of 
flow in an exterior domain found in |15| . It was shown there that, if div u £ (t, •) = 0, 
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curl u £ (t, •) = u) e (t, •), and if u £ (t, •) is tangent to r e and vanishes at infinity, then 
there exists unique a = a(t) £ M. such that one can write u £ (t, •) as: 

u E (t,-) = K £ [uj £ (t,-)} + a{t)H £ ; 

(see Section 3.1 of JH] for details). For the initial data 2JI we have, of course, 
a(0) = a. We will prove that a(t) — a for any t > 0. This fact relies on a result 
whose prove we defer to Section 8, see Corollary^] Using the notation introduced 
above we have: 

u £ - o% = K £ [uj £ ] + a(t)[l - <p £ ' X °}H £ + [o(t) - a}<p £ > x °H £ . 

It will be proved in Corollary [TBI that u £ — o§ belongs to Lf£ c ({0, oo); L 2 (U e j), 
although the estimate blows up as s — > 0. With this result in mind we first note 
that, from the estimates for K £ proved in ^S], we expect to have that G 
Lfo c ([0, oo); L 2 (n e )). In fact, if oj £ were compactly supported in space for each 
fixed time then this would follow from estimate (2.8) in |15| : the adaptation to our 
case is possible but escapes the scope of the present work. Next we recall that the 
harmonic vector field H £ is smooth (because the conformal map T £ is smooth and 
extends smoothly to T £ ) and has 0(1/ \x\) behavior near infinity. Therefore we find 
that [l-ip e ' X °]H e e L 2 (n e ), but tp £ ' X °H £ <£ L 2 (Tl e ). Hence the only way for u £ - o% 
to be square-integrable is for a(t) — a = 0, as we wished. 

Since the flow u £ satisfies the no-slip condition at any positive time, the circu- 
lation around r e at t > vanishes. We make use once more of Stokes' Theorem to 
conclude that = a(t) — m £ (t), where m £ (t) — J n o-i £ (t, x) dx. We can now account 
precisely for the mass of vorticity produced at the boundary in the initial layer. We 
have: 

m £ (t) = a. 

4. Initial data asymptotics 

The purpose of this section is to study the limit, as e — > 0, of the initial velocity 
fields Uq. We begin by introducing some notation. 

For each function / defined on II e , we introduce Ef, the extension of / to K 2 , 
by setting Ef = in efl. 

Lemma 8. If f G Wj ' c (II e ) and if its trace vanishes on the boundary T £ then 
E f G W ?oc and E ^f = VEf. 

The proof of this fact is elementary and we leave it to the reader. 
We will now introduce notation which will be used in the remainder of this paper. 
We denote by P the Leray projector on all of K 2 . Additionally, we introduce the 
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cutoff 

(4.1) rf = rf{x) = ^ s {x) = <p{\T £ {x)\) = <p(\T(x/e)\), 

where T £ , ip £ ' £ and ip, were introduced in Section 4, see (|3.1|) . I|3.4fl , Note that there 
exists a constant C > independent of e such that r/ £ (x) = 1 in {|x| > Ce}. 
Let 

( 4 - 2 ) = ^xf 

be the kernel of the Biot-Savart law in all of R 2 , / i— > JC* /. Note that we denoted 
the same vector field by H in 12.41) . The different notations used for the same vector 
field are natural since x 1 - / (27r|a;| 2 ) plays two very different roles - one as the kernel 
for the Biot-Savart law for the full plane and another as a harmonic generator for 
the cohomology of the exterior of any disk centered at the origin. 
Let Uq be as in (|3.2f) and 7 as in (|3.5() . 

Lemma 9. Let uq = K, * luq + . Then we have that 

¥[ifEu £ ] -> u m T>'(M. 2 ) 

as e — * 0. 

Proof. We split u £ in a different way than before: 

u £ Q = (K £ [lu ] + mH £ ) + 1 H £ = v £ + jH £ , 

where m was defined in l|3.5() . By Lemma 4.2 in ^S], r/ £ EH £ — > strongly in 
Lj 1 oc (IR 2 ) as e — > 0, and by Lemma 4.1, in rj £ EH £ is divergence-free, so that 
Pt] £ EH £ = i] £ EH £ . Therefore, 

7 P[?f -> 7i? in X>'(R 2 ). 

All that remains to prove is that V[t] £ Evq] ^ K, * ujq in 2?'(]R 2 ). To see this, we 
begin by observing that t) £ Evq is uniformly bounded in L°°(R 2 ), see Theorem 4.1 in 
|15j . Furthermore, we have additional control over the behavior of Evf, at infinity, 
so that there exists a constant C > 0, independent of e, such that \ti £ Evq\ < C/\x\. 
This follows from the explicit expressions for K £ , H £ given in (3.5) and (3.6) of 
|15j . from estimate (2.8) in |15| and from the compactness of the support of ojq. 
Therefore, t] £ Evq is also uniformly bounded in L P (M. 2 ) for all p > 2. Fix p > 2 and 
let C e L^QR 2 ) n LP(R 2 ) be a weak-* limit of {r] e Ev^}. 

Next we observe that div t) £ Evq — Wt) £ -Evq and curl t] £ Evq = V^-rf -EvQ+rfujQ. 
The cutoff if is such that \Vr] £ \ is bounded by C/e and supported on a set of 
measure Ce 2 . Thus, Wrj £ (and V^rf) converges to zero strongly in L q (U. 2 ) for any 
1 < q < 2. Hence, div rj £ Ev^ — > and curl t] £ Evq — > ojq strongly in L 1 . This, 
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together with the convergence of a subsequence of rfEv^ to £ weak in LP implies 
that div ( = and curl £ — uq in the sense of distributions. Using that £ £ LP for 
p < oo, we obtain that £ = /C * c^o- 

Since we identified the limit, we have actually proved that rfEv^ — 1 /C * loq 
weakly in L p , without the need to pass to subsequences. Therefore, as P is linear 
and continuous from LP to itself, it follows that 

F[t) e Evq] -± F[K *lj }=1C*ujo, 

which concludes the proof. 

□ 

5. The impulsively stopped rotating cylinder 

In this section we will study the small obstacle asymptotics in the special case 
where the obstacle is a small disk and the initial flow is harmonic. The purpose 
of this discussion is just to motivate and illustrate our main result, as this special 
case is not rigorously required for the remainder of the analysis. 

Let us begin with a physical interpretation of our problem. Consider an infinite 
solid cylinder of radius r > immersed in a viscous fluid occupying the whole 
space outside the cylinder. If the cylinder rotates with constant angular velocity 
A, boundary friction will induce rotational motion in the surrounding fluid which, 
one expects, will settle to a steady flow with velocity uq of the form 

Xr 2 x x 
U0 = ^' 

see PP for a discussion of this example. 

We consider viscous flow in the exterior of the cylinder with initial velocity 
uq, imposing the standard no-slip condition u = at |z| = r for positive time. 
Physically, this corresponds to first "preparing" the initial data by rotating the 
cylinder for a long time, letting the flow settle into the steady configuration uq, 
and then suddenly halting the motion of the cylinder at time t = 0. A shorthand 
description of this situation is that of the flow induced by an impulsively stopped 
rotating cylinder. The inconsistency between initial and boundary data in this 
problem generates a rather singular initial layer in the fluid motion. 

This problem has both translational symmetry, along the axis of the cylinder, 
and rotational symmetry about the same axis. If we assume that these symmetries 
are preserved by the flow then the translational symmetry allows us to reduce the 
problem to two dimensions, and the rotational symmetry cancels the nonlinearity 
and further reduces the dimension, so that the equations of motion reduce to the 
Stokes equation in the exterior of a disk. One may find in pQ an explicit treatment 
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of this problem, involving passing to polar coordinates, using separation of variables 
and expressing the solution by means of Fourier-Bessel integrals. 

The problem we wish to address next is the small obstacle limit of the impulsively 
stopped rotating cylinder as posed above. This means that we consider II £ = {|x| > 
e}. We have H £ — x- 1 / (2ir\x\ 2 ), independent of e. In the notation of the previous 
section we pick ujq — and 



(5.1) 



X 

Note that, 7 = 2tt\e 2 , so that fixing the circulation 7 independent of e means that 
the angular velocity A of the obstacle must blow up as the obstacle becomes smaller. 

We consider u £ = u £ (x,t) and p £ — p £ (x,t) solving Ij3.3|l with initial data as 
above. It is a nice exercise, which we leave to the reader, to prove that the solution 
preserves circular symmetry. Preserving the symmetry means that the velocity 
remains tangent to concentric circles about the origin, with the pressure and the 
magnitude of velocity both invariant under rotation. One consequence of circular 
symmetry is that u £ ■ Vu £ is a gradient field, so that we can absorb the nonlinearity 
into the pressure term. Therefore, u £ satisfies the Stokes system on II £ . 

We introduce the Lamb-Oseen vortex as the unique solution U of 

U t + U ■ VU = -Vp + uAU, in (0, 00) x R 2 
div U = 0, in [0, 00) x M 2 

lim^oo U(t, x) = 0, for all f > 

see ^| ^] . We are now ready to state and prove the main result in this section. 

Theorem 10. Let 7 G M. Then the extension of velocity, Eu £ , converges weakly 
m L 2 OC ((0,^),L 2 OC (R 2 )) to jU. 

We will present an outline of a proof for this result, highlighting the main ideas, 
for the sake of illustration. In Section 10 we will give a complete proof of our main 
theorem, which includes this example as a special case. 

Proof. By linearity we can assume, without loss of generality, that 7=1. 

Let us begin by collecting the estimates for Eu £ . As in Lemma|51 we decompose 

u e = b e + i s + o e = + (1 - ip)H £ + ipH £ , 

with ip the same cutoff used in (|3.4fl . 
We have: 

(1) Eu £ is bounded in L^ o ° c ((0, 00); LP(R 2 )) and in L^ oc ([0, 00); LP(R 2 )), for any 
2 < p < 00, q = 2p/ (p — 2) and 1 < q < q ; 
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(2) Eu e is bounded in L^ c ((0, oo); L 2 oc (R 2 )); 

(3) VEu £ is bounded in L 2 OC ((0, oo); L 2 {R 2 )) and in Lf oc ([0, oo); L 2 (M 2 )), for 
1 <p< 2; 

(4) £w e - <^iJ is bounded in L£ c ((0, oo); L 2 (R 2 )). 

Indeed, we can use Lemma to get the first item. The second and third items 
come from Lemma |5] Theorem ^ and Proposition [2J The last item comes from 
Lemma and Proposition 

Let uj e = curl u e be the vorticity. Let tp = ip(x) G C^°(R 2 ) and 9 = 9(t) G 
C£°((0, oo)) be test functions and, for each <5 > e consider </j 5 = ^(|a:|/#), with y> 
as before. The vorticity satisfies the vorticity equation in Ii E , which reduces to the 
heat equation in this special case, by symmetry. Multiplying the vorticity equation 
by ip s ip9 and integrating by parts we find 

(5.2) f [ tp s iP0 t Euj E dxdt = -v [ [ 6Eu e A((p s ip) dxdt. 

Jo Jr 2 Jo Jr 2 

Fix Euj £k a weakly converging subsequence in L 2 OC ((0, oo); L 2 (R 2 )) and let u> be 
its weak limit. Let u be the weak limit of Eu £k , so that u> = curl u for t > 0. Note 
that all the a priori estimates (1-4) above pass to u and u. One may pass to the 
limit as £k — > in identity (|5.2|l to obtain 



(5.3) / / tp s 6 t -4>uj dxdt = -v \ I 9cjA(ip s iP)dxdt. 

Jo Jr 2 Jo Jr 2 

Next, we assume that "0(0) = 0. 
We have 

A(<p 5 ip) = i/jA<p s + 2VV- • V(f S + <p s A<ij>. 

Note that ip s — 1 — * in LP for any 1 < p < oo, and that both V<// and ipAip 5 
converge weakly to zero in L 2 as 6 — > 0. The last convergence follows from the 
hypothesis ip(0) = 0. Therefore, taking the limit S — > in (|5.3I) we obtain 



/>OC /> /-OO P 

(5.4) / / 6 t ipuj dxdt = -v I I 9uiAipdxdt. 

Jo Jr 2 Jo Jr 2 

This is enough to prove that the associated velocity u is a weak solution of 

the Navier-Stokes system in the full plane. To see that, take <£> a divergence-free, 

compactly supported test vector field in the plane and define ip = A -1 curl This 

means that <£> = V^ip. For each R > let x R = (1 — vd^l/-^)) an d define 

$ B = V ± ((^-^(0)) X i? ). 
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Take i\i — i(jr — (ip — ip(0))x R m <|5.4J) . We integrate by parts to obtain 

/>oo /> />oo /> 

(5.5) 0=/ / 6 t <5> R udxdt- v\ \ duAip R dxdt = h - vl 2 . 
Jo Jr 2 Jo Jr 2 

We wish to pass to the limit R — ► oo in both Ii and I 2 ■ First note that 



h = / / M(Ax K )W> - ^(°)) + 2V * • v ^ + X A$ cfecft. 
Jo Jr 2 

Observe that both Ax fl and V% fl converge weakly to in L 2 and that \ R 
pointwise. Therefore, 

lim I2 = I / BluA-i/} dxdt. 
R ^°° Jo Jr 2 

On the other hand, since ipH is independent of time, 

h= I ! (G t $R)(u-(p(\x\)H)dxdt 

JO JR 2 



" 



Therefore, 



pOC p poo p 

lim h = / / [0 t $){u - <p{\x\)H) dxdt = / / 6 t <S>udxdt 
Jo Jr 2 Jo jr 2 



Putting together the limits of I\ and I2, 15.5(1 and integrating by parts in the 
vorticity term we obtain 



t <&udxdt = -v I I 9A$udxdt. 
10 JR 2 Jo Jr 2 

This is the weak formulation of the velocity equation, under the symmetry as- 
sumption. As for the initial data, it is not difficult to prove that the initial data for 
the limit problem is the limit (in the sense of distributions) of the initial data for 
the approximating problems. This can be shown either by taking the test function 
9 not vanishing in and including the initial data in the weak fomulation and in 
the process of passing to the limit, or showing that the velocity converges uniformly 
in time (up to time t — 0) with values in Hr^ in the same way as in Proposition 
im Therefore, by Lemma the velocity u satisfies u(0,x) = H. From uniqueness 
of weak solutions, see ^JEI, the conclusion follows. 

□ 

A result analogous to Theorem EH holds for solutions of the Stokes problem with 
arbitrary initial data, not necessarily symmetric, i.e., for any given 7, the small 
obstacle limit of solutions of the Stokes problem with initial data of the form (|3.2|) 
is the solution of the Stokes problem in the full plane with initial data given by 
Lemma El The key issue here is linearity. As we will see in the next section, the 
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smallness condition on the initial circulation is needed to deduce estimates in the 
initial layer, but it is solely due to the nonlinearity. 

Finally, we would like to use the impulsively stopped rotating cylinder as an 
illustration of what is taking place in the initial layer. We describe the events 
at time t = for the flow generated by an impulsively stopped rotating cylinder 
noting that, for an (infinitcsimally) small positive time, the fluid velocity vanishes 
at the boundary but has a nonvanishing limit as one approaches the boundary from 
inside the fluid. Tangential discontinuities in fluid velocity are called vortex sheets 
in hydrodynamics. The effect of impulsively stopping the rotation of the boundary 
amounts to placing a vortex sheet at the boundary and letting it diffuse into the 
bulk of the fluid through viscosity. The same rough picture describes what happens 
in the initial layer for our general problem as well. 

6. INITIAL-LAYER AND THE NONLINEAR EVOLUTION 

We have fixed an arbitrary initial vorticity wo and circulation 7, and hence we 
must understand solutions of the Navier-Stokes equations with initial data which 
does not satisfy the no-slip boundary condition. The effect of the consequent initial 
layer can be understood roughly as that of placing a vortex sheet at the boundary 
T e and letting it evolve, diffusing into the flow. In Section 5 we saw that, if we 
consider linear evolution as described by the time-dependent Stokes semigroup, 
the a priori estimates from Section 2 are enough to establish the small obstacle 
asymptotics. However, the problem of resolving this initial layer for the full Navier- 
Stokes system and obtaining uniform estimates for the small obstacle problem is 
rather more delicate and it is the subject of the present section. 

Let u £ be the solution of the Navier-Stokes equations with initial velocity 
u% given by (O and l et v e = K] as in Section 3. Let W s = u £ ~ v £ . Let P £ 
be the Leray projector on LT e . The evolution of W £ is described by the following 
system: 

(6.1) W £ - i/P £ AW e + ¥ £ div (W £ <g> W £ + W £ ® v £ + v £ ® W £ + v £ <g> v £ ) = 0, 

with the initial condition VF e (0, x) = and the boundary condition W £ = on T e . 

We introduce the weighted-in-time norms. Let p > 1 and / : (0,T) — > L P (H £ ) 
measurable. Let T > 0. We use the following notation: 

\\f\\ P ,T= sup i*~*||/(V)IU*(rie)- 

0<t<T 

The use of these norms for the Navier-Stokes equations was pioneered by H. Fujita 
and T. Kato, see for example 
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Lemma 11. Let 2 < p < oo. There exist positive constants Co and C p such that, 
ifO<T< C v 3 and \y\ < Cqv then 

p + 2 

\\W e \\^ T <C p v— , 

for every < e < Eq. 

Proof. We use Duhamel's principle to write 

W E (t) = - [ Sl{t - t)P £ div (W £ <g) W £ + W £ <g> v £ + v £ <g> W £ + v £ ® v £ )(t) dr. 
Jo 

Take the L p -norm and apply Theorem^ estimate (S3), with 2 < q < p to obtain 

H^lli^n.) < / \\S s u {t-T)¥ e div (Ty e ®^ e +^ e ^ £ +w £ 0VK £ +w e ^ £ )(r)|| iP( n e )dr 
Jo 

i-t 



<K 3 / M*-T))"' + ? " 5 (l|W £ (T)IU*i||W re (T)IU« +l|W re (T)|| iW ||» e (T)||£« + 

Jo 

+ \\v £ (T)\\ L « 1 \\v e (T)\\ Lq2 )dT, 

where q\ and qi are chosen so that 1/q = l/qi + l/<?2 & n d we have used Holder's 
inequality. Next we use the definition of the (p, i)-norm to find 

\\W £ \\ LP{ n c) < ^ 3 (||^ £ || gi , t ||M/ £ || g2 , t + ||^ £ || 9l , t |K|| 92 , t + 

+ \\v e \\ qi ,t\\v% 2 ,t) f{v{t-r))-^-\A- l dr. 
Jo 

We note that, for any a > —1, (3 > —1, we have 
(6.2) f {t- T ) a T f3 dr <C{a,l3)t a+fs+1 . 



Jo 

The proof of this inequality is an elementary calculation. 

We wish to use l|6.2[) with a = —1/2 + 1/p — l/q and = 1/q— 1. Assume that: 

(6.3) 2<q<p. 

Note that this condition implies 

a= -- + ---> -1 and /3 = - - 1 > -1. 

2 p q q 

Therefore, we find 
(6.4) 

\\w%, t < Cv-^-H\\w £ \\ qi , t \\w £ \\ q2 , t + \\w% ltt \\v% 2 , t + \\v% ut \\v% 2tt ) 

We divide the remainder of the proof in two steps: p — 4 and any p > 2. 



2D INCOMPRESSIBLE VISCOUS FLOW AROUND A SMALL OBSTACLE 



19 



First assume p = 4. Set qi = qi = 4, so that q = 2 and (|6.3(1 is satisfied. In this 
situation, Hfi.4[l gives that X(t) = \\W £ \\4 t t satisfies 

(6-5) * 2 +(lKlk 4 -^)*+IKHl t >0. 

Note that X(0) = and the parabola described by i)6. 5fl has, at t = 0, two 
distinct nonnegative roots. We observe that, as long as this parabola has two distinct 
nonnegative roots r\(t) < ^(t), we have that inequality l|6.5[l together with the 
continuity in time of X , n and T2 imply that 

(6.6) < X(t) < n(t). 

The condition for the polynomial above to have two distinct roots is 



(6-7) M| w e| kt __j _4 M 2 it>0 . 

Since ||w e ||4,t > and ^ 3 / 4 /C > we find that I6.7|l is equivalent to 
(6-8) Klkt < "3^- 

Furthermore, under the above assumption, the two distinct roots are also nonneg- 
ative. The (4, i)-norm is nondecreasing in t and hence, in order to guarantee that 
the polynomial in l|6.5|) have two distinct nonnegative roots, it is enough to verify 
(|6.8I) for t = T. Now we use the linear estimates from the previous section to find 
conditions under which (|6.8() is valid at t = T. 

First, recall that v £ = b e + o £ + i e , and that both bg and Oq belong to L P (I1 £ ), 
for p > 2, with L p -norms uniformly bounded in e, see Lemma We use estimate 
(SI) from Theoremnwith p = q = 4, together with Lcmma[S]to deduce: 

IKIU,t< sup ti\\b e (t)+o £ (t)\\ LHnE) + sup ti\\i e (t)\\ LHne) 

0<t<T 0<t<T 

< tfiTiaibSIU^.) + ||og|U4 ( n £ )) + K\-y\v-l < C(Ti + \ 7 \u~i). 
Choose Co > so that the conditions 

(6.9) T < C is 3 and \j\ < C v 

imply with t = T. 

Assuming now that (E3 are valid, using l|6.6|l we have that 

(6.10) \\W e \U, t = X(t) < ri (t)<Cv 3 /\ 

for < t < T. This concludes the proof in the case p = 4. 

For any p > 2 we bootstrap the (4, T)-estimate in the following way. We return 
to Hfi.4|l and set qi = q 2 = 4. We then impose (j6.9|l to obtain 
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||W e || p , T < Cv-i + $-H\\W S \\l T + ||M/ e || 4 ,T||f e ||4,T + \\V £ \\1 T ) < C&)v*&. 

□ 

We conclude this section with the observation that 

||W" e ||2,T = ||W -e ||x,«.(( 0l T);i^(ne))) 

so Lemma ITTI actually provided a renormalized energy estimate on the initial layer. 

7. Global-in-time nonlinear evolution 

In the previous section we obtained a priori estimates for W £ in the initial layer 
which are uniform in e. We will now splice the information we already possess with 
a standard energy estimate, in order to obtain a result which is global in time. We 
retain the context introduced in the previous section. 

Lemma 12. Letl<p<2. Then W £ G L% c ([0, oo); L 2 (n E ))flL 2 oc ((0, oo); H 1 ^)), 
W e G Lf^QOjOo);.// 1 ^)), and the respect ive norms are bounded independently of 
e. 

Remark 13. Note that the bound in Lf oc ((0, oo); H 1 ^)) means that W £ is 
bounded in L 2 ((S, T); H 1 (II e )) for any < 5 < T, but not necessarily for 5 = 0. 

Proof. We rewrite the evolution equation (|6.1|l for W £ as 

W£ - vAW £ + (W E + v E ) ■ VW £ + W e ■ Vv e + v e ■ Vv E = -Vp, in (0, oo) x IT e 
div W £ = in [0, oo) x II e 
r(0,-)=0on{f = 0}xn E 
W e {t, ■) = on [0,oo) x T e 

We multiply the equation above by W e and integrate to obtain 

S = ~\\W e f L * + v\\VW e \\ 2 L2 = - [ [W £ ■ (W £ ■ Vv e ) + W e ■ (v £ ■ Vv £ )} dx 
2 dt ' J n _ 

= f [v £ ■ {W £ ■ \7W £ ) + v £ ■ (v £ ■ \7W £ )} dx 

< ll^ £ |U4||VM/ £ || L 2||^|| L 4 + \\VW £ \\ L 2\\V £ \\ 2 L ,. 

We will use the following interpolation inequality: 

\\w £ \\ L , < c\\w £ \\% 2 \\vw s \\% 2 , 

with a constant C > independent of e. This inequality in the case of R 2 can 
be found in Chapter 1 of |18| . To obtain the corresponding inequality in LT £ , one 
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simply extends W e to K 2 by setting it identically equal to zero inside eSl. As W £ 
vanishes on r e , the extension has i/^-norm in the plane identical to the i/^-norm 
of W £ in il e . Finally one uses the inequality in M 2 on the extension. 
We proceed with the estimate of £ : 

£ < c\\w^ 2 2 ^w s \\% 2 \\v e \\ L , + nvwi^KH^ 
< |liv^ni 2 + ^\\w s \\Uv £ \\i* + l\\v £ \\U 

where we used Young's inequality to estimate each of the products above. Next, 
we use Lemma [S] to deduce 

C 



kill. < u1 



Hence, 



±\\w-\\h+^\\h<^w-\\h + £ t , 



for some constant C independent of e. Gronwall's inequality now gives, for any 
< h < t 2 , 

W £ (t2,-)\\h , f t2 \\VW%s,-)\\ 2 L2 ^ v 2 v 2 \\W%h,-)\\ 2 L , 



C 7 ' 1 ) t c/^ +V J t s c l vi ds ~ t c/v 4 t c/u4 + t c/vl 

First choose t\ = Co^ 3 /2, with Co given in Lemma HT1 It follows from Lemma ITT1 
with p = 2 that 

ll^ e (ii,-)l!i, <cv 2 . 

Therefore, 

) ' 



(7.2) \\W £ {t,-)\\ 2 L2 <Cv 2 



ct\ l 



for any t>t\, and we conclude that W £ is uniformly bounded in L?° ([0, oo); L 2 (H £ )) 
as desired. 



Next, we return to l|7.1|l . for the derivative estimate. Let a > 0, multiply (|7.1|) 
by i a ^ c l v ) 1 arLC j integrate the resulting inequality with respect to t\ from to 
<2- We obtain, 

(7.3) t s a \\VW s (s, -)f L2 ds < a ^4^ \-t a 2 + T S*- 1 \\W%s, OH!, da . 

JQ v . ° JO 



Since we already know that W e is uniformly bounded in Lf£ c ([0, oo); L 2 ), this esti- 
mate implies that W s is bounded in £ 2 OC ((0, oo); H 1 ), uniformly in e. Moreover, if 
1 < p < 2, then the choice a = (2 — p)/2p above allows to conclude that W £ is also 
bounded in Lf oc ([0, oo); H 1 ). □ 
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8. Velocity estimates 

In this section we derive global estimates on velocity using the analysis performed 
thus far. Before we begin, we require the following interpolation inequality. 

Lemma 14. Let 2 < p < oo, qo = 2p/(p ~ 2), 1 < q < qo. Let r > q(p — 2)/p and 
set 

rp - q(p - 2) ' 

If p > 2q assume further that r < q{p — 2)/(p — 2q). Then 9 > 1 and for any 
interval ICR and any f G L r (I; H 1 ^ 2 )) H L e (I; L 2 (R 2 )), we have 

\\f\\Li(I;LP(M^)) < C 'll/llL'-(7 ; ) ^(R2))ll/lli' / /(/ ;L 2( R 2 ) )- 

Proof. We start by recalling the following standard interpolation inequality: for 
any g £ H 1 (R 2 ) we have 

(8.1) ||fl|Uvci-)(R») < C\\g\\ H , < qi^ll^'Hflll^, for any < s < 1. 

Fix exponents p, q, r and 9 as in the statement of this lemma. Observe that, 
if p > 2q then 9 > 1 if and only if r < q(p — 2)/(p — 2q); we hence assume this 
further restriction on r if p > 2q. In the other case, p < 2q, there is no additional 
restriction on r to guarantee that 9 > 1. 

Next, hx an interval ICR and let / € L r (I; H 1 (R 2 )) n L e (I; L 2 (R 2 )). Let 
s = (p — 2)/p, so that 2/(1 — s) = p. We use (|8.1ll and Holder's inequality to obtain: 

ll/ll W ;LP(R2)) = I \\f(r, OHL < C jf ||/(r, -)ll^r 2)/?, |l/(T, -)f// P <*r 

<- /-ill f II 2 )/P II f ||2g/p 

— ^lu \\L r (I;H 1 (S. 2 ))WJ Hl 8 (7;L 2 (R 2 ))' 

which concludes the proof. The condition r > q(p — 2)/p was used in Holder's in- 
equality when estimating the product of two functions in I/pMp- 2 ) (j) an( j j^vl (2<?) ( 
above, so as to guarantee that rp/q(p — 2) > 1. 

□ 

Theorem 15. Let it &e the solution of <|3.3fl initial velocity Uq as in (|3.2|) 

and recall that Uq = + ig + Oq. TTien £/ie following hold true. 

(1) Let 2 < p < oo, go = 2p/(p — 2) and 1 < q < q . Then {Eu e } is bounded 
m L* c ((0,oo);L p ) n L q loc ([0,oo)\L p ). 

(2) TTie family {Eu £ — Eoq} is bounded in Lf£ c ((0, oo); L 2 ) and the family 
{Eu e - Eo £ ~ Ei £ } is bounded in L% c ([0, oo); L 2 ). 

(3) The family {Eu e } is bounded in L^ c ((0, oo); Lf oc ). 

(4) For any 1 < p < 2, we have {V Eu £ } is bounded in Lf oc ((0, oo); L 2 ) n 
LL([0,oo);L 2 ). 
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Proof. Statement (1) involves two estimates: the first one on the open time interval 
(0,oo) and the second on the closed interval [0, 00). We begin by addressing the 
first estimate. 

Fix 2 < p < 00. Fix < 6 < T and set I = (5,T). We first show that 
{Eu E - Eo £ } is bounded in L°°((5, T);L 2 (R 2 )) n L 2 ((S, T); H 1 (R 2 )). We write 

(8.2) Eu s - Eo e Q = E(i e + h e ) + E{o e - o e ) + EW £ = A t + A 2 + A 3 . 

We observe that A\ is bounded in L°°(I; H 1 (R 2 )). To see that, choose 1 < r < 2 
and use Theorem ^ together with Lemma Eland Remark to obtain 

(g 3) (l|i e ll^(n e) + l|b £ || L2 (n e )) + t* (||Vi E || L 2 (rie) + || Vb £ || i2(rie )) 

< (ifi + ^2)(||i§|Ur(n.) + ||b§|U, (n .)) <K{r), 

for some K(r) > 0, independent of e. The estimate on A\ follows from Lemma |S1 
together with the inequality above. For A 2 , we use Proposition |3 together with 
Lemma [S] to conclude that 

(8.4) \\o e (t,-)-o e \\l 2(Jle) + v f ||Vo £ ( s ,-)-Vo ||| 2(rii) d s < ^t||Vo ||| 2(ne) < C, 

Jo 

for some C > independent of e. This, together with Lemma |S1 implies that A 2 
is uniformly bounded in L°°(/; L 2 (R 2 )) n L 2 (I; H 1 ^ 2 )). For the estimate on A a , 
we simply use Lemma 1121 together with Lemma |H1 showing that A3 is uniformly 
bounded in L°°(I;L 2 (R 2 )) n L 2 (I: H^R 2 )) as well. 

We use Lemma IT41 with q — qo and r = 2, so that 9 = 00, to conclude that 
{Eu E - Eo E } is bounded in L qo (I; L P (R 2 )). Next we note that {Eo E } is uniformly 
bounded in L P (M. 2 ) by Lemma \5\ which concludes this portion of the proof. 

We now address the second part of (1), which is an estimate on the closed time 
interval [0, 00). The difficulty here is that we do not have Leray-type estimates on 
the pieces of u s all the way down to t — 0, so that the result becomes more delicate. 

Fix 2 < p < 00 and 1 < q < qo. Let T > and set I — [0, T]. We consider again 
the decomposition Ij8.2|l and we estimate each piece. Estimate l|8.H|l . together with 
LemmaElimplies that A 1 is uniformly bounded on L ri (I; H^R 2 ))!!^ 1 (I; L 2 (R 2 )), 
for any 1 < r\ < 2 and any 1 < Q\ < 00. We use Lemma Il4l with p and q as above. 
We need to find re [1,2) satisfying the restrictions in LemmaHHin order to be able 
to use r = r±. This is always possible because the restriction on r always includes 
t > <?(p — 2)/p, and q < q is equivalent to q{p — 2)/p < 2. This implies that A\ is 
bounded in L q (I; L P (R 2 )). For A 2 , we merely observe that (|8.4|l gives an uniform 
bound in L°°{I; L 2 (R 2 )) n L 2 (I; H 1 ^ 2 )) on A 2 , which in turn yields the desired 
estimate. To treat ^3, we put together Lemma ITSl and Lemma |H1 to conclude that 
A 3 is uniformly bounded in L°°(/; L 2 (R 2 )) n L r2 (I; H\R 2 )), for any 1 < r 2 < 2. 



24 



IFTIMIE, LOPES FILHO AND NUSSENZVEIG LOPES 



Clearly, this is enough to obtain the estimate in L q (I; L P (M. 2 )) for A3. The proof 
of (1) is concluded once we recall the observation that {Eoq} is uniformly bounded 
in L P (M. 2 ), which we already used in the proof of the first part of (1). 

We now address statement (2), which also consists of two estimates. The proof 
of the first estimate in (2) is contained in the proof of the first part of (1), the 
estimate on the open time interval. As for the second estimate in item (2), we 
write 

Eu £ - Eo% - Ei £ = Eb £ + E{o £ - of,) + EW £ . 

We have already shown that the second and third terms in the decomposition above 
are bounded in Lf£ c ([Q, 00); L 2 (M 2 )). The first term satisfies 

\\Eb s (t, -)\\ L , < \\b £ (t, Oilmen.) < C\\bl\\ LHUe) < C, 

by Theorem ^and Lemma 

The third item, statement (3), can be obtained from (2) by observing that, by 
Lemma[S]and Lemma|S| Eoq is uniformly bounded in L r (R 2 ), for any r > 2, which 
is contained in L 2 oc (R 2 ). 

Statement (4) again consists of two estimates, one on the open time interval, the 
other on the closed interval. The estimates on the open time interval are trivially 
contained in the proof of the first estimate in item (1), once we observe that Lemma 
0and Lemma|Hlgive a uniform estimate in L 2 (R 2 ) for V-Eoq. Similarly, the proof 
of the second part of (4) is contained in the proof of the second part of (1), together 
with the L 2 (K 2 ) estimate for V-Eoq which we have just derived. 

This concludes the proof. 

□ 

The last result in this section is an estimate for "fixed e", which was already 
used to deduce that the amount of vorticity generated at the boundary in the 
initial layer, m £ (t) — J n w £ (t, x) dx is the same for any e and equals a. 

Corollary 16. For eachfixedO < e < e we have: Eu £ ~Eo £ Q £ L^ c ([0, 00); L 2 (U e )). 

Proof. This is an immediate consequence of item (2) in Theorem 1 151 together with 
the observation that E\ £ £ Lf£ c ([0, 00); L 2 (H e )) for each fixed e. This last fact 
follows from the estimates on the Stokes semigroup in Theorem ^ and the fact that 
EXq £ L 2 (Jl s ) with an L 2 -norm that blows up as | loge|i, see Lemma |5] □ 

9. Compactness in space-time 

As far as a priori estimates go, the last ingredient we require is uniform control 
on how solutions evolve in time. This is often very easy to accomplish once spatial 
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estimates are in place because ultimately the PDE itself is nothing more than an 
expression of time-derivatives of the solution in terms of spatial information. In 
our case, however, there are two difficulties that will make this step of the analysis 
somewhat involved: (i) the spatial estimates available are for Eu £ , which does not 
satisfy a PDE: and (ii) the nonlinearity in our problem is quadratic, but the estimate 
up to time zero in Theorem 1151 item (1) is L p , p > 2, which entails problems at 
infinity. We deal with these difficulties through the following main ideas: we use the 
vorticity equation to describe the time evolution, we use the interplay of vorticity 
and velocity and we renormalize problem terms. 

Let $ be a smooth, compactly supported vector field and consider P the Leray 
projector for the plane. We consider the Hodge decomposition of the vector field 
<&, given by $ = P<3> + ($ — P$). The divergence- free part P<I> is smooth, but not 
compactly supported. In fact, 



see the proof of Proposition 1.16 in [22]. Let ip = ip(x) be the stream function 
associated with P$, so that V tp — We assume that ip(0) = 0, at the expense 
of having ijj = 0(1) at oo. Clearly, |^(a;)| < \x\ || V-0||l°° , so that using the Sobolev 
imbedding H 2 L°° followed by the fact that Vip = V(A)" 1 V ± • $, a zeroth 
order singular integral operator acting on $, we obtain 



We now observe that for 1 < q < oo and \ € C^°(M 2 ), there exists a constant 
C = C(q, x) > such that 



Indeed, as P is a zero-th order singular integral operator we have 

||P[x$]|U»(fi2) < C||x$Hl*(r2) < C , ||x||l 5 (r2)||$||l-(r2 ) < C(g,x)||$||ij2(R2)- 

Recall the extension operator E, introduced in beginning of Section 4. For each 
e > 0, consider the cutoff rf introduced in (|4.1|) . We are ready to state and prove 
the main result in this section. 

Proposition 17. The sequence {P[r] e Eu E ]} is precompact in Lf£ c ([0, oo);H^(R 2 )). 
Proof. Fix <!> a smooth, compactly supported vector field and let 



(9.1) 




(9.2) 



(9.3) 



l|]P[X $ ]IUo(K 2 ) < C||*||i?=(R=)- 



^ = ^( x ) = [(A)" 1 curl $](x) - [(A)" 1 curl $](0), 
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satisfying (|9.2|) . For each t > 0, we introduce an auxiliary functional F £ — F e (t) E 
H- 2 (R. 2 ) denned by 

(F e (t), $) = / (Eu £ (t 7 x) - Eo%{x)) ■ (V-V )(a;)V>(20 dx. 



The proof will be divided into two steps. We will show that, for each (tx,t2) C 
[0,oo), we have {F[rj £ Eu £ } + F £ } is bounded and equicontinuous as a function of 
(ti,t 2 ) intoiJ ; ~ 2 (R 2 ) and we will show that F £ -> strongly in L% c ([0, 00); tf- 2 (R 2 )). 
The desired conclusion follows from these two steps by using Arzela-Ascoli's The- 
orem. 

Let us begin by proving that F £ — > strongly in Lf£ c ([0, 00); 7J~ 2 (R 2 )). Indeed, 
we use Theorem (|9.2|) and the properties of the cutoff rf to deduce 

(9.4) \(F e (t),$)\ < \\Eu s (t,-) - Eo e \\ L 4iPV x r) s \\ L , 
<C\\Eu £ (t r )-Eo £ Q \\ L J f N 2 ||$|!^ 2 |V^| 2 ^ 

\J\x\<Ce J 

<Ce\\Eu £ (t r )-Eo £ a \\ L 2\m H2 

< Ce(\\Ei £ (t, OH^+H^Ct, ■)-Eo e -Ei%t, -)\\ L 2)\m H 2 < Ce(C 1 \\oge\^+C 2 )\\n„2, 

by Lemma |3 Theorem ^ and Theorem 1151 Clearly, this proves our assertion. The 
proof of the other assertion is a bit more involved. 

We introduce a cutoff for infinity. For each R > 0, let x R — X R { X ) — ^~ ( p{\ x \/R)- 
The vorticity lo £ = curl u £ satisfies the equation 

oj £ t + u £ ■ Wuj £ = isAui £ . 

Let < <i < t 2 < 00 and denote the interval [ti, t 2 ] by J. We multiply the vorticity 
equation by rj £ ipx R i integrate in space and time between t\ and t 2 , and integrate 
by parts to obtain 

(9.5) J [Eu £ (t 2 , ■) - Eu e (h, •)] ■ V x (rf V>X fl ) dx 

= / {Eu £ ■ VEuo £ )-q £ ^x R dxdt -v \ \ {AEu; £ )r] £ ^x R dxdt = h~h- 
Jt x J Jt 1 J 

We first estimate I\ . We integrate by parts and deduce 

J (Eu £ ■ Vri £ )Eu £ ^x R dxdt + J J {Eu £ ■ (F^) A -)Eui £ ri s x R dxdt 

- J J J {Eu £ • Vx fl )£wW dxdt = -I n + I u - J 13 . 
Using Holder's inequality first in space and then in time we have 

I/11I < f 2 ll^ll^n^ll^llW IMIx fl IU~ dt 

Ju 
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< C||VVr ? E || L 4|| J Bu £ || i 3 (J . i 4 ) ||i;a; e || i9 /5 (j ;. i2) |t2 - *i| 1/9 < C||*||fli|*2 - t^ 9 , 
where in the last inequality we used Theorem 1151 and we used again 19.2J1 along 
with the properties of rf . 

Similarly, we have 

|/ia| < IU*l|S« e |U»(j i i*)ll^w e || z ./. (Ji£a) |ta ~ *i| 1/9 < C\\n H ,\t 2 - t^l\ 
since P is a zeroth order operator and H 2 L 4 . Finally, 

|A 3 | < CH^IU- \\Vx R \\L4Eu s \\ L3iJ;Li) \\Euj s \\ l2/3(j . l2) < Cim^R- 1 / 2 , 

as 

|V(oo)| = II(A)- 1 curl $](0)| <C\\<P\\ H 2 

and S7x R — 0(1/ R), supported on a set of measure 0(R 2 ). 
Therefore, 

(9.6) limsuplAl < C\m H2 \t2-h\^ 9 . 

R->oo 

Next we treat I 2 . We integrate by parts and use the fact that the support of 
Vrj e and of V\ R are disjoint to obtain: 

h = 



Euj £ (Ar] £ )ijx R dxdt + v / / Euj e r) e (Aip)x R dxdt 
Jti J 

Eu) e rfi)(&x R ) dxdt - 2v J J Euj £ (x R Vr] e + ? ? £ Vx K )(P$)- L dxdt 
= vl 2 \ + vl 22 + vl 2 3 - 2vl 2i . 
By arguments similar to those used for I\ we have 

ij 21 i < \\&^nL4E^\\ L ^ [ j., L i ) \t 2 - tii 4/9 < c\m H *\h - ^i 4/9 ; 

|/ 22 | < ||A^|U= / 2 \\Euf\\& dt < C\\ curl *|| i3 |i a - ^l 4 / 9 < CH^IIflalta - ti| 4/9 ; 
Jt 1 

|/ 23 | < \\Ax r \\ L 2U\\l™ f 2 \\Eu; s \\ L z dt < C\\$\\ H 2; 
Jti 

< WmL~\\x R V V E + V e Vx R \\l* / 2 ||^ £ || L 2*<C||<i>||^|i2-ii| 4/9 . 

Jt! 

Therefore, 

(9.7) limsup|I 2 | <C||$|| H2 |*2-ti| 4/9 - 

R — >oo 



We expand the left hand side of identity l|9.5|) to find 

[Eu £ {t 2 , •) - Eu £ (t u •)] • V ± (V £ ^X R ) dx 
[Eu E (t 2 , -)-Eu e (h, -)KV VWX* dx+ f[Eu £ (t 2 , -)-Eu e (h, -)W(V^) X R dx 
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+ J [Eu s (t 2 ,-)-Eu s (h,-)}-r } e ij{V x x R )dx = A 1 +A 2 + A3. 

We will show that each of the Ai's has a limit when R — > oo. To see that, first 
note that 

Eu £ (t 2 ,-)- Eu%t!,-) 

= [Eu £ (t 2 , •) - Eo £ - Ei £ (t 2 , ■ )] - (i i , • ) - Eol - Ei £ (t i , • )] + Ei £ {t 2 , • ) - Ei £ (f i , •) , 

which belongs to L 2 (M. 2 ), for each fixed e > and < t\ < t 2 < oo. To see this 
note that the first two terms are bounded in L 2 by Theorem ^] whereas the last 
two terms were estimated in L 2 , with a logarithmically growing norm as e — > 0, in 

llIEIIl- 
Therefore, since V -1 ^ and V ± ip = P<& are both square integrable functions, it 
follows by the Dominated Convergence Theorem that 

lim Ax= I [Eu £ (t 2 , •) - Eu £ (t u •)] • (V ± if)ipdx 

R^oo J 

and 

lim A 2 = [[Eu s (t 2 ,-)-Eu e (t l ,-)]-r ] s {V X 'il))dx. 

R-*oo J 

Furthermore, it is easy to see that V^x^ converges to zero weakly in L 2 when 
R — > oo. As [Eu £ (t 2 , •) — Eu £ (ti, •)] • ifip does not depend on R and belongs to L 2 , 
we infer that A3 — > as i? — > 00. 

We have found that the left hand side of identity 1)9. 5|) has a limit as R — > 00. 
We can rewrite this limit as follows 

lim (Ai + A 2 + A 3 ) = (P[?f £u £ ](i 2 , •) + F £ (i 2 ), $) - (Pfo e £!u e ](ti, •) + F e (i x ), $). 

— >oo 

On the other hand, by identity Q9.5|l . and using l|9.6[) and l|9.7[) we have 
(9.8) lim \A 1 +A 2 + A 3 \<\m-isup\I 1 \ + \I 2 \<C\\<S>\\ H 2\t 2 -t 1 \ 1 / 9 , 

which shows that ¥[ifEu £ ] + F £ is equicontinuous as a function of time into H~ 2 . 

We conclude this proof by showing that {P[q £ Eu £ ] + F £ } is uniformly bounded 
in L°°(J; ff^fK 2 )). We do not need to prove that {F £ } is bounded in this space 
because we have already shown that F £ — > as e — > in L°°(J; if ~ 2 (K 2 )). The only 
thing left is to prove the boundedness of {P[?f To this end, let x € C^°(R 2 ), 

fixp> 2, 1 < r < 2 and write 

\( X V[ifEu £ },<P)\ = \( V £ Eu £ ,¥[ X m 

< \(rf(Eu* -Eol - Ei*),F\x*})\ + \(rf Eo* ,F\x*])\ + \(r, e Ei £ , P[ X $])| 
< H^^-— ^og— ||lP[^<l>]|| i2 + ||^og|| ZjP ||]P[^cl>]||^ / ^_ 1) ^||^i-|| i ^||]P[^<l>]|| i ^ /( ,_ 1> 

<C(x,P,r, J)\\*\\h*, 
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where in the last inequality we used Theorem^ LemmaEl Theorem l 1 51 and relation 
(|9.3(l . Note that C is independent of t E J. 

It follows from Arzela-Ascoli that, for each [ii,t2] C [0, oo) and each Br c R 2 , 
there is a subsequence of ¥[rj £ Eu £ ] which converges strongly in L°°([ti, £2]; H~ 3 (Br)) 
By taking diagonal subsequences we may assume that there is a subsequence which 
converges strongly in £j£ c ([0, 00); H^(R 2 )). 

This concludes the proof. 

□ 

Remark 18. It follows from the proof of Proposition El that any strong limit of 
{P[ V £ Eu £ }} in L~ ([0,oo);i^ c 3 (R 2 )) in fact belongs to C([0, 00); i^ c 3 (R 2 )). This 
is true as {¥[r] e Eu £ ] + F e } is equicontinuous and bounded as a function of time 
into -ff;~ 2 and therefore, by Arzela-Ascoli, its limits are continuous. Furthermore, 
F e 0, so that the limits of {¥[n e Eu 6 ]} and of {¥[n E Eu e ] + F E } are the same. 

10. Passing to the limit 

In this section we state and prove our main result. Let us begin with an improve- 
ment of the space-time compactness we have, which is a consequence of Proposition 
1171 obtained by means of interpolation. 

Lemma 19. The sequence {Eu 5 } is precompact in L 2 oc ((0,oo) x R 2 ). 

Proof. By Lemma|Sl first order derivatives of functions that vanish on T £ commute 
with the extension operator, and therefore, for any positive time, Eu e = P[Eu £ ]. 
We write 

Eu e = P[(l - tj^Eu 5 } + P[n e Eu £ ] =B X + B 2 . 

First we note that B\ — ► strongly in Lf oc ([0, 00) x R 2 ). Indeed, let us fix < t\ < 
t-2 < 00 and set J = [ti,t2]. By Theorem properties of the cutoff rf and the 
fact that the Leray projector is continuous from L 2 to itself we have 

||P[(1 - 7f )£ W e ]|| i2(JxR2) < C||(l - r, £ )^ £ || i2(JxR2) 

< C||l - ^|| L 4 (R 2)||£;u £ || L 2 (J . L 4 (R2)) < Ce 1/2 , 

which proves the desired estimate on B\. 

Next we work on B 2 . We know from Proposition El that B2 is precompact in 
L% c ((0,oo);Htol(M?)). We will show that, for any 1 < q < 2, B 2 is bounded in 
L 2 loc ((Q,oo); Wi*(R 2 )). The result will follow by interpolation. Fix 1 < q < 2 and 
let q* = 2q/{2-q) > 2. By TheoremlTBI {Eu e } is bounded in lf oc ((0, 00); L 4 (M. 2 )). 
Since |?7 e | < 1 and since F is continuous from L 4 (R 2 ) into itself, it follows that 
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B 2 is bounded in Lf oc ((0, 00); L 4 (R 2 )) which can be continuously imbedded into 

LL((0,oo);LL(K 2 ))- 

What remains is to show that derivatives of B 2 are also uniformly bounded 

in L 2 OC ((0, 00); L' oc (R 2 )). Since the gradient and the Leray projector P are both 

Fourier multipliers, the gradient commutes with P. Therefore, 

D{P[ifEu £ }) = ¥[ri s (DEu e )} + ¥[(D v £ )Eu £ ] = B 21 + B 22 . 

By Theorem H5l DEu £ is bounded in Lf oc ((Q, 00); L 2 (R 2 )). Since |?7 e | < 1 and 
P is continuous from L 2 (R 2 ) to itself, we immediately obtain the desired esti- 
mate for B 2 \. As for the term B 22 , we use Theorem 1151 once again to obtain 
that Eu e is bounded in L 2 fJ {q '~ 2) ((0, 00); L«* (M 2 )) and we recall that Drf is 
uniformly bounded in L 2 (M. 2 ). With this, we have that (Drj £ )Eu £ is bounded in 
L^ C V(? *" 2) ((0, 00); L«(R 2 )), continuously imbedded into i 2 oc ((0, 00); L 9 (R 2 )). This 
concludes the proof. 

□ 

We will prove that limits of the sequence {Eu £ } are solutions of the Navier- 
Stokes equations in a suitable weak sense. To be precise, we formulate the notion 
of weak solution we will use. 

Definition 20. Let u G L 2 OC ((0, 00) x R 2 ) n C([0, oo);V'(R 2 )). We say that u is a 
weak solution of the incompressible Navier-Stokes equations with initial velocity uq 
if, for any divergence-free test vector field \& £ C£°((0, 00) x IR 2 ), we have 

poo f 

A*) dxdt = 0. 

Jo Jr 2 

Furthermore, for every t > 0, div u(t, •) = in the sense of distributions and 
u(t, •) — uq in the sense of distributions as t — > + . 

Recall that JC denotes the kernel of the Biot-Savart law, as introduced in . 
We are finally ready to state and prove the main result of this work. 

Theorem 21. Any strong limit u of {Eu £ } in L 2 OC ((0, 00) x R 2 ) is a weak solution 
of the incompressible Navier-Stokes equations in R 2 with initial velocity given by 

Uq = K * LUq + jH. 

Remark 22. As {Eu £ } is precompact, by virtue of Lemma ll9l there exists at least 
one such strong limit. 

Proof. For each e sufficiently small, choose < S < 1 such that {\x\ > 26} C n e . 
Clearly, if {|a:| > 25} C IT S0 then {\x\ > 26} C II £ , for all e < e . Also consider 



2D INCOMPRESSIBLE VISCOUS FLOW AROUND A SMALL OBSTACLE 



31 



R > 2 > 25. We use the cutoff ip introduced in Section 4 to define: 

ip s = /(z) = p(|a;|/a) and X R = X R (x) = 1 - <KM/fl). 

As in the proof of Proposition El we let $ be a smooth, compactly supported 
vector field in R 2 , which, in addition, we assume to be divergence- free. We define 
t/> = ip( x ) = [(A)" 1 curl - [(A)" 1 curl $](0). Recall that V 1 ^) = $ and that 

■0 satisfies @£2J. We also consider 6> = 0(t) e C~((0,oo)). 

We use the test function ip s 0ipx R i which belongs to C£°((0, oo) x II e ) in the weak 
form of the vorticity equation. We can rewrite the integrals on II £ as full plane 
integrals using the extension operator to obtain the following integral identity 

poo p poo P 

/ / Euj £ 9 t <f 5 ^x R dxdt+ / / Eu e Euj e ■ dW(ip 5 ^x R )dxdt 
Jo Jr 2 Jo Jr 2 

(10.1) +v f ( ELu £ 6A(ip 5 iJjx R )dxdt = 0. 

Jo Jr 2 

Our first step is to pass to the limit e — > in this identity, while keeping S and R 
fixed. Let u be a strong limit in Lf oc ((0, oo) xt 2 ) of a subsequence Eu £k of Eu £ . Wc 
observe that Ew £k — > curl u = u strongly in L 2 OC ((0, oo); 7J ; ~^(R 2 )). Similarly, we 
may also deduce that u is divergence- free in the sense of distributions. The passage 
to the limit is immediate in the linear terms of HlO.lf) . For the nonlinear term we 
recall that, by Theorem IT51 {Ew £k } is uniformly bounded in Lf oc ((0, oo); Z 2 (R 2 )). 
Hence a subsequence of {Euj £k } converges weakly in L 2 OC ((0, oo); L 2 (R 2 )). Using 
the convergence Etu £k — > ui strong in £ 2 OC ((0, oo); ^/^(R 2 )) and uniqueness of weak 
limits we conclude that Ecu £k — w weakly in £ 2 OC ((0, oo); L 2 (R 2 )), without passing 
to further subsequences. Now Eu £k Eu £k is a weak-strong pair, so that we can pass 
to the limit in the nonlinear term as well. We arrive at the identity 

pOO /> />oo /> 

Ji + J2 + J 3 =/ / uj8 t v s i/jx R dxdt+ / uuj ■ 0X7(tp s ifjx R )dxdt 
Jo Jr 2 Jo Jr 2 

(10.2) +u f [ u6A( i p s ipx R )dxdt = 0. 

Jo Jr 2 

Now we pass to the limit both 6 — > and R — > oo in each separate term in 
(EE2JI. We begin with J x . 
First we observe that 

(10.3) u-acp H^u-Fe £&(((), oo); L 2 (R 2 )), 

where a was introduced in l|3.2|) . (3 in Lemma [3] and Ao in Lemma and H is the 
harmonic vector field introduced in l|2.4(l . Identity (|10.3(l follows from the conver- 
gence Eu £k — > u, from the fact that i?w e - Eo £ Q is bounded in Lf% c ((0, oo); L 2 (R 2 )), 
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see Theorem 1151 and from the fact that Eoq — > atp ^4!jr) H strongly in L 2 (M. 2 ) by 
Lemma [S] 

Next, write u> = V • u and integrate by parts to obtain 

poo p poo P 

Jx = - / u-9 t \' ± {<p s ipx R )dxdt = - / (u-F) ■ t W- L (ip s il;x R ) dxdt, 
Jo Jr 2 Jo Jr 2 

where we have used the fact that F does not depend on time, so that the additional 
integral vanishes. 
We write 

poo p poo P 

Ji = - / (u — F) ■ 9 t (W ± ip s )^ X R dxdt -I (u-F)- 9 t tp s <S>x R dxdt 



- f [ (u-F)- t <p s il>(y ± X R ) dxdt = -J u - J 12 - J 13 . 
Jo Jr 2 

It is easy to see that Vip s converges to zero weakly in L 2 (R 2 ) when 5 — *■ and 
therefore, 

(10.4) lim lim J u = 0. 

R^oo S^O 

On the other hand, Vx R also converges to zero weakly in L 2 (R 2 ) when R — > oo. 
Furthermore, ip s — > 1 poinwise when 8 — > 0, so by dominated convergence, we find 
that 

(10.5) lim lim J 13 = 0. 

H->oo <5^0 

We also note that <p s x R converges poinwise to 1 as 5 — > and R — > oo (no matter 
which order), so that, by dominated convergence, we deduce that 

poo P P OO P 

(10.6) lim limJi2= / / (u — F) ■ 9 t <& dxdt = / / u -9^ dxdt. 

R^COS^O J J R2 J J R 2 

Putting together ltTO|) . fRFty and fSJHfy we obtain 

(10.7) lim lim Ji = - / / u -9^ dxdt. 

R^oc 5^0 J J M 2 

Next we treat the nonlinear term J2. First note that the uniform estimates on 
Eu e contained in Theorem 1151 imply that u € Lf oc ((0, 00); L 4 (R 2 )). The argument 
is the same we used to prove oj e L 2 OC ((0, 00); L 2 (R 2 )) 

We write 

pOO P pOO P 

J 2 = / uuj -9(Wip s )^x R dxdt- / uv •6<p s $ ± x R dxdt 
Jo Jr 2 Jo Jr 2 



23- 



+ / f uu • 9(p s ip(V X R ) dxdt ee J 2 i - J 22 + J; 
Jo Jr 2 

We have that 

POO _ 

\J*i\< / \6\\\u\\ L 4u\\ L 4^<p 5 \\ L *dt = 0(V6) 
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and 

\Jaa\ < Ml* 
We conclude that 



pOO 

/ \6\\\u\\ L 4u\\ L 4V X R h*dt = 0(R- 1 / 2 ). 
Jo 



(10.8) lim lim(J 21 + J 23 ) = 0. 

R^oo S-^0 

In addition, by dominated convergence we have that 

p oo p roc p 

(10.9) lim IimJ 22 = / / uuu ■ 6<S>^ dxdt = - I I (u ■ V)u0 ■ $ dxdt, 
fl->oo<5->-0 Jo Jr 2 Jo Jr 2 

where this last equality follows from the identity u ■ Vw — (uui) = V(|it| 2 /2), 
together with the fact that <£> is divergence free. 

Therefore, using 1)10.8)1 and (|10.9(l and integrating by parts we find 

(10.10) lim lim J 2 = - / / [(u ■ V)6>$] • udxdt. 

fl->oo 8^0 J J R 2 

Lastly we treat J 3 . Once again, we write 

/>OC /> />00 f- 

J 3 = v / / uj6(A<p 5 )<tpx R dxdt + v / u 6(p 5 (A^) X R dxdt 



ui6ip d ip(Ax M ) dxdt + 2u / ujOV^ ■ {(Vtp d )x U + V (^X )) dxdt 
Jo Jr 2 

= J3I + J32 + J33 + ^34- 

Using, similarly to what we have already done, that: u> € Lf oc ((0, 00); L 2 (R 2 )), 
ipA(p s weakly in L 2 (R 2 ) as 5 -> 0, Ax R -> strongly in L 2 {R 2 ) as R -> 00, 
Vip s weakly in L 2 (R 2 ) as <5 -> 0, Vx fl weakly in L 2 (R 2 ) as R -> 00 and 
dominated convergence, we deduce that 

(10.11) lim lim(J 31 + J 33 + J 34 )=0. 
Therefore we obtain, integrating by parts, 

pOO p poo p 

(10.12) lim \imJ 3 = is / w 0(Aip) dxdt = -v \ \ u-6A$dxdt, 

R.^OoS^O Jo J R 2 J Q J R 2 

Recall that J x + J 2 + J 3 = 0, so that, adding flu"?)) with 1)10. 10|l and l|10.12[l we 
find 



(10.13) 



/ / (u-6 t $ + [(u-V)6$\-u + ini-9A$)dxdt = 0. 
Jo Jr 2 



We observe that linear combinations of products of smooth, compactly supported 
functions of the form 9& are dense in C£°((0, 00) x M 2 ). With this observation 
and (|10.13fl we find that u satisfies the integral identity in Definition [201 We have 
already noted that u is divergence- free in the sense of distributions. All that remains 
is to show that u S C([0, 00); V'(R 2 )) and that u(t, •) u in V as t -> 0. 
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First note that in the proof of Lemma IT§I we showed that Eu £ — P[q £ Eu £ ] — > 
strongly in L? oc ((Q, oo) x R 2 ). Therefore, W[rf»Eu ek ] -> u in Lf oc ((0, oo) x R 2 ). By 
Remark 1181 we have a subsequence of {F[rj £k Eu £k ]} which converges strongly in 
£^ c ([0,oo);ff~ c 3 (R 2 )) to a limit v G C([0, oo); ff^R 2 )). It follows by uniqueness 
of limits (in L 2 OC ((0, oo); HZ^), for example) that u(t, ■) — v(t,-) for almost all 
t £ (0, oo), which implies that u can be identified with v. This in turn implies that 
u € C([0, oo); 2?'(R 2 )). Furthermore, as F[rj £k Eu £k ) converges to u uniformly in 
time with values in HT^t, one has that F[r] £k Eu £ k ] converges to uq in H7^. On the 
other hand, LemmaElsays that F[rj £k Eu £ k ] converges to K, *ojq +"fH in the sense of 
distributions. By uniqueness of the limit in T>' , we conclude that uq = K, *ljq + jH. 
This concludes this proof. □ 

Remark 23. At the end of the proof above we showed that the initial data for the 
limit problem is attained in C([0, oo); 2?'(R 2 )). We can actually prove a stronger 
statement, namely that there exists a positive constant C > such that 

\(u(t)-u ,<t>)\ < Cll^lff^ 1 / 9 , 

for all <j> e V. Indeed, to see this 6x<f> € C C °°(R 2 ) and let ip = tp(x) = [A" 1 curl <j>](x)- 
[A" 1 curl 0](O). Consider the sequence of approximations {u £k } constructed in the 
proof above. For each Ek recall the auxiliary functional F £k used in the proof of 
Proposition El given by 

(F £k (t),(j)) = J {Eu £k (t,x) - Eo £ k (xj) ■ (V A -rj £h )(x)i>(x)dx. 

Write 

(10.14) (u(t) - u (h & = L! + L 2 + L 3 + L 4 , 

where 

L 1 = (u(t) -F[if k Eu £k }(t) - F e *(t),<j>); 

L 2 = (F[7f k Eu £k ](t)+F £k (t)-F[r] £k Eu £ k } - F £k (0), (/>}; 

L 3 = (F £k (0) ,<!>); 

L 4 = {F[r) e *Eu e k ]-u ,<f>). 

By Rcmark llSI we have that lim efc ^o L\=Q uniformly in time. By the argument 
in the proof of Proposition^! see estimate (|9.8|) . we find that \L 2 \ < CH^H/pt 1 / 9 . 
By estimate ()9.4|l we know that \\m ek ^L 3 = 0. Finally, by Lemma El we obtain 
that lim E(c ^o L4 = 0. Therefore, using l|lU.14fl we deduce 

\(u(t)-u ,<fr)\ < CIMI^t 1 / 9 , 

as desired. 
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Remark 24. The solution obtained at the limit is much more regular then what 
is required for a weak solution as stated in Definition 1201 For example, the a priori 



estimates given in Theorem 1151 imply immediately that u 6 ((0, oo); L p ) for all 
p E (2,oo) and Vu € Lf oc ({0, oo); L 2 ). 

Our result above provides strong compactness of viscous flows around a small 
obstacle but does not address actual convergence. The passage from compactness to 
convergence is clearly reduced to the issue of uniqueness for the limit problem. The 
issue of uniqueness of solutions for the 2D incompressible Navier-Stokes equations 
with initial data of the form uq is classical and delicate. Let us briefly review the 
related literature. The first relevant results are due to G. Benfatto, R. Esposito 
and M. Pulvirenti, see who showed uniqueness for initial flows of the form 
J2liH(x — Xi) if 5D|7»| is sufficiently small, and to G.-H. Cottet, see |1|, who 
showed uniqueness for small initial vorticities which are general bounded measures. 
Later, Y. Giga, T. Miyakawa and H. Osada generalized this uniqueness result for 
initial flows of the form K,*lj with u) a Radon measure with sufficiently small atomic 
part, see [141 116j . As we observed in the Introduction, this smallness condition is 
closely related, in a technical sense, to the smallness condition on 7 which we also 
had to impose, see (|6.9JI . Recently the uniqueness assumption on the atomic part of 
oj was removed, first by T. Gallay and C. E. Wayne for initial flow of the form jH 
in |12| . see also |10j . and then for general K. * loq initial flows with cog an arbitrary 
Radon measure by I. Gallagher and T. Gallay in UJ. These results are a byproduct 
of the remarkable large-time asymptotics results obtained by Gallay and Wayne in 



From the point of view of the present work, the natural question is whether Gal- 
lagher and Gallay's result in |Uj implies uniqueness of the weak solutions obtained 
as small obstacle limits in Theorem |21 Gallagher and Gallay's proved that there 
is at most one solution v — v(x,t) to the 2D Navier-Stokes equations with initial 
vorticity a Radon measure satisfying 

(1) the vorticity w = curl v E C((0, T);L 1 (R 2 ) n i°°(IR 2 )). 

(2) |m-,£||z,i < K for all t > 0. 

(3) u>(-, t) — ujq + 760 as t — * 

(4) the vorticity w is a mild solution of the vorticity formulation of Navier- 
Stokes equations on any time interval [to, t\] compactly contained in (0, 00), 
i.e., 



EH 




V-e' 




for any < to < t < t\. 
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In our result, we obtain a weak solution u in the sense of Definition 1 with 
initial vorticity loq + jS. The additional regularity stated in Remark 1241 should 
be more than enough to prove that the solution u is a strong solution for t > 0, 
thus placing our result in the setting of Gallagher and Gallay. For example, the 
condition u S £* c ((0, oo); L P (M. 2 )) places our solution on the Serrin criterion curve, 
see |25| . which is a standard condition for uniqueness of weak solutions, see also 

m 

In short, there is no uniqueness result stated in the literature which includes 
precisely our solution, but clearly such a result is expected and should be easy to 
prove. We will not address this issue any further, as it escapes the main purpose 
of this paper. 

11. Conclusions 

The purpose of this section is to interpret what we have done in a broader context 
and to point out some directions for improvement and further work. Our basic 
problem was to find conditions under which the presence of a single small obstacle 
could be ignored in the modelling of large scale flow. The precise formulation we 
used, working in the unbounded exterior domain and fixing the large scale flow 
by choosing an initial vorticity ujq and a circulation 7, was convenient from the 
mathematical point of view, but it was far from physically natural. 

There are many ways in which to formulate mathematically the problem of 
placing a small obstacle in a given incompressible viscous flow. To be more precise, 
let us suppose that we are given a smooth background flow Ub — Ub(t,x) and say 
we wish to insert a small circular obstacle (centered at the origin) within this flow. 
That would mean adding a correction u £ c = u £ c (t, x) so that the new flow u £ = Ub+u £ c 
satisfies the no-slip boundary condition at \x\ — e. It makes sense to assume, for 
simplicity, that there are no more boundaries in the problem and that Ub vanishes 
at infinity. We also expect u £ c to be sharply localized, and hence it should vanish 
at infinity as well. In this context, there is no reason for the background vorticity 
ujb = curl Ub to have vanishing total integral. Let us denote by uib the total integral 
of tub- The vanishing of mb is equivalent to requiring the background flow to have 
globally finite kinetic energy, a physically reasonable assumption. However, infinite 
energy flows associated with smooth but mass-unbalanced vorticity have often been 
considered (see Section 3.1.3 in [22] for a discussion), so we do not assume m& = 0, 
but we consider this case to be of special physical relevance. 

Let cu £ = curl u £ . As we have seen, there exists a — a(e) such that u £ = 
K £ [uo £ ] + a(e)H, where H was defined in 12.4(1 . If u £ should converge to Ub then 
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j oj € — a(e) — ► 0. If this were not the case, then the kinetic energy of u £ would 
blow up near the origin, a complicated consequence of the method of images for the 
disk, see Moreover, since both — m^H] — 0(l/\x\ 2 ) and \u e — a(e)H\ = 

0(l/\x\ 2 ) at infinity (see [22]), we must have a(e) — » mb as e — ► for the kinetic 
energy of the difference Ub — u e to be bounded at infinity. 

In other words, at any fixed time, if a flow u e around the obstacle approximates 
the background flow Ub without infinite energy discrepancies either at infinity or 
locally, then it must have vorticity ui £ with total integral close to nib and its har- 
monic part should also be close to mt>. It is therefore very natural to model this 
problem with a solution of the Navier-Stokes equations for which: 

(1) the initial vorticity loq has integral approaching m&, 

(2) ojq approximates u)b(0,x) in L \ 

(3) the initial circulation 7(e) = a(e) — J luq — > as e — > 0. 

The main point of this discussion is to claim that the results obtained in our article 
can be easily adapted to include physically natural approximations such as those 
described above. Therefore, although we chose a very specific way of formulating 
the small obstacle limit, reasonable alternatives should lead to the same result. 

From the discussion above and the physical interpretation of the case 7(e) ^ 0, 
given in Section 5, together with the natural scaling of this problem, we can see 
that the cases 7(e) = or 7(e) — > are by far the most interesting situations. 
However, there would have been no substantial simplification of the argument by 
restricting our problem to the case 7(e) = oil). Moreover, the situation in which 
7(e) = O(l) is mathematically very interesting. Indeed, the smallncss condition on 
the initial circulation only appears when we assume 7(e) = 0(1). Also, there is a 
discrepancy between the results obtained for the inviscid and viscous cases when 
7(e) = O(l), which suggests that the limits e — ► and v — > do not commute. 

We have assumed throughout that ujq was smooth and compactly supported. 
How much regularity on loq did we really use? The answer is none. We actually 
needed uo bounded in Lf oc and L 2 '°° and nothing else. We contrast this with the 
inviscid argument, where we needed loq in L p , p > 2. 

Let us turn to some problems which arise naturally from our work. One par- 
ticularly interesting question is the issue of considering both the viscosity and the 
obstacle small. This should be a difficult problem, because the wake due to an 
obstacle becomes more pronounced and turbulent as viscosity vanishes. It is well- 
known that, for full plane flow, one can take the vanishing viscosity limit, obtaining 
solutions of the incompressible 2D Euler equations, see OEll^- In the presence of 
a material boundary, the vanishing viscosity limit is a classical open problem, even 
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if the flow is very smooth. The difficulty is due to the boundary layer. The prob- 
lem of taking the vanishing viscosity limit outside a very small obstacle interpolates 
nicely between the full plane result and the open problem of taking the vanishing 
viscosity limit in the presence of a fixed material boundary. In fact, this question 
is one of the main motivations of the present work and it is still under consider- 
ation by the authors. Taking into account the result obtained in this paper it is 
clear that one should first pursue the small viscosity problem in the case 7(e) = 
or 7(e) = o(l) as e — > 0, since the smallness condition in our convergence result 
gets more restrictive as viscosity vanishes. With this future work in mind we have 
included the specific dependence of our estimates on viscosity for as long as it was 
practical. 

A second problem is to extend our analysis to velocity fields which are constant 
at infinity, in order to include the classical case of a material body moving in a fluid 
with roughly constant speed. 

Yet another problem that arises from our work is to remove the smallness con- 
dition on the initial circulation. The parallel between our convergence problem 
and uniqueness for the limit flow suggests a strategy. Is it possible to adapt the 
entropy-entropy flux techniques used by Gallay and Gallagher for the uniqueness 
problem to the small obstacle asymptotics? 

A fourth problem is to obtain an asymptotic description of the correction term in 
the small obstacle limit, i.e. a description of the "wake" associated with the small 
obstacle. Finally, one can consider a whole host of related problems, described 
loosely as the study of limit flows in singularly perturbed domains. For instance, 
one can study limit flows in a bounded domain with one or more small obstacles, 
or in a domain composed of a small neck joining two fat domains, or in a domain 
having a long thin tail, etc. 
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